Dynamics of domain walls in magnetic nanostrips 
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We express dynamics of domain walls in ferromagnetic nanowires in terms of collective coordi- 
nates generalizing Thiele's steady-state results. For weak external perturbations the dynamics is 
dominated by a few soft modes. The general approach is illustrated on the example of a vortex 
wall relevant to recent experiments with flat nanowires. A two-mode approximation gives a quanti- 
tatively accurate description of both the steady viscous motion of the wall in weak magnetic fields 
and its oscillatory behavior in moderately high fields above the Walker breakdown. 



Dynamics of domain walls in nanosized magnetic wires, 
strips, rings etc. is a subject of practical importance 
and fundamental interest 0, Q • Nanomagnets typically 
have two ground states related to each other by the sym- 
metry of time reversal and thus can serve as a memory 
bit. Switching between these states proceeds via cre- 
ation, propagation, and annihilation of domain walls with 
nontrivial internal structure and dynamics. Although 
domain-wall (DW) motion in macroscopic magnets has 
been studied for a long time Q, new phenomena arise 
on the submicron scale where the local (exchange) and 
long-range (dipolar) forces are of comparable strengths 
0| . In this regime, domain walls are textures with a rich 
internal structure @, As a result, they have easily 
excitable internal degrees of freedom. Providing a de- 
scription of the domain- wall motion in a nanostrip under 
an applied magnetic field is the main subject of this pa- 
per. We specialize to the experimentally relevant case of 
thin strips with a thickness-to- width ratio t/w <C 1. 

The dynamics of magnetization is described by the 
Landau-Lifshitz- Gilbert (LLG) equation [|| 



m = 7H c ff x m + a m x rh. 



(1) 



Here m = M/|M|, H off (r) = -SU/5M(r) is an effec- 
tive magnetic field derived from the free-energy func- 
tional [/[M(r)], 7 = g\e\/2mc is the gyromagnetic ratio, 
and a < 1 is Gilbert's damping constant Q. Equa- 
tion |T]) can be solved exactly only in a few simple cases. 
Walker S] considered a one-dimensional domain wall 
m = m(x, t) in a uniform external magnetic field H||x. 
At a low applied field the wall exhibits steady motion, 
m = m(x — vt), with the velocity v ~ "fHA/a, where 
A is the wall width. Above a critical field H w — aM/2 
magnetization starts to precess, the wall motion acquires 
an oscillatory component and the average speed of the 
wall drops sharply. Qualitatively similar behavior has 
been observed in magnetic nanostrips however, nu- 
merical studies demonstrate that Walker's theory fails to 
provide a quantitative account of both the steady and 
oscillatory regimes Q. 



We formulate the dynamics of a magnetic texture in 
terms of collective coordinates £(t) = {£o, £i, • • •}, so that 
m(r, t) = m(r, {£(<)})• Although a magnetization field 
has infinitely many modes, its long-time dynamics — most 
relevant to the motion of domain walls — is dominated by 
a small subset of soft modes with long relaxation times. 
Focusing on soft modes and ignoring hard ones reduces 
complex field equations of magnetization dynamics to a 
much simpler problem. In Walker's problem, the soft 
modes are the location of the domain wall and the pre- 
cession angle; the width of the wall is a hard mode [2, @] . 
Partition of modes into soft and hard depends on char- 
acteristic time scales, determined e.g. by the strength of 
the driving field. 

Equations of motion for generalized coordinates {£(£)} 
describing a magnetic texture can be derived directly 
from the LLG equation {!]). They read 



Gijij + Fi - \ ,,£, - 0. 



(2) 



Here = —dU/d^i is the generalized conservative 



force conjugate to £j, while Ti 



and G%j 



-Ga are 



the damping and gyrotropic tensors with matrix elements 
described below. The three terms in Eq. ((2|) can be traced 
directly to the three terms in the LLG equation ([1]). 

To derive Eq. © , take the cross product of Eq. (Q]) with 
m and express the time derivative of the magnetization 
in terms of generalized velocities, m(r, £) = (3m/9^)^, 
to obtain 
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(3) 



Here J = hqM/j is the density of angular momentum. 
Taking the scalar product with dm/d^i and integrating 
over the volume of the magnet yields Eq. $2$ with 

= - f SU/Sm-dm/d&dV = -dU/d&, 
ry(£) = a,J J dm/dti-dm/dtjdV, 
Gijit) = J J m • (dm/d^ x dm/d£j) dV. (4) 
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FIG. 1: Top: A model of the vortex domain wall proposed 
in Ref. [Tol . Dashed lines denote Neel walls emanating from 
the topological edge defects. Bottom: Absorption and re- 
emission of the vortex at the edge. Note the reversal of the 
polarization p of the vortex core. 



Eqs. and (01 generalize Thiele's result [9j for steady 
translational motion of a texture to the case of arbitrary 
motion. 

We apply this general approach to the dynamics of the 
vortex domain wall [Bj], a texture that consists of three 
elementary topological defects: a vortex in the bulk and 
two antihalfvortices confined to the edges 111] . A strong 
shape anisotropy forces the magnetization into the plane 
of the strip, with the exception of the vortex core 
Soft modes of the wall are associated with the motion of 
these defects, and we start with a model parameter- 
ized by the (X, Y) coordinates of the vortex (Fig. 1). In 
low applied fields, the wall exhibits translational motion 
that can be described by a single collective coordinate 
£o = X, representing the softest (in fact, zero) mode 
with an infinite relaxation time tq = oo. At higher driv- 
ing fields the steady motion breaks down and the vortex 
core exhibits oscillations in both longitudinal and trans- 
verse directions accompanied by slow drift along the strip 
2]. An additional dynamical variable £i = Y, is required 
to describe the dynamics. The new mode has a finite 
relaxation time t\. In the vortex domain wall the char- 
acteristic time of the motion is time T it takes the vortex 
to cross the strip. When 



n>_>I}, >T> T 2 >T l > 

soft 



(5) 



hard 



the soft modes £o and £i must be treated as dynamical 
variables. All other modes are hard; they adjust adia- 
batically to their equilibrium values. As the driving field 
increases, the vortex moves faster and eventually T will 
become shorter than the relaxation time of the next 
mode, at which point the two-mode model will break 
down. While tq is infinite due to translational symmetry 
of the wire, t\ is also long because of the special kinemat- 
ics of vortex cores (see discussion below). If t\ ^> T2 we 
expect to have a substantial range of driving fields where 
the two-mode approximation applies. 



Next we discuss the general aspects of the dynamics 
in the one and two-mode regimes. We approximate the 
potential energy U(X,Y) by its Taylor expansion to the 
second order in X and Y: 



U{X, Y) w -QHX - xrQHY + kY 2 /2 



(6) 



The X dependence comes in the form of the universal 
Zeeman term —QHX, where Q — 2fioMtw is the mag- 
netic charge of the domain wall independent of the ex- 
act shape of the texture. Zeeman force also pushes the 
vortex in the transverse direction, which is reflected in 
the linear in Y term, dependent on the vortex chirality 
X = — 1(+1) for clockwise (counterclockwise) circulation. 
This term is consistent with the lack of y i— > —y reflection 
symmetry; the numerical coefficient is r « 2. The trans- 
verse restoring potential kY 2 /2 comes from the dipolar 
and exchange energies. 

The antisymmetric gyrotropic tensor Gxy = —Gyx = 
4irqJt reflects a special topology of the vortex core, 
namely its nonzero skyrmion charge [l3| 



q = (1/4tt) J m ■ (d x m x d y m) d 2 r 



np/2, (7) 



where n = +1 is the 0(2) winding number and p — 
M z 1 1 M z | — ± 1 is the out-of-plane polarization of the core 
[lj] . A vortex core moving at the velocity V experiences 
a gyrotropic force F 9 = pGz x V, where G — 2nJt is 
the gyrotropic constant. The equations of motion ([2]) for 
two dynamic modes read 



^xx ^xy - pG 
Tx y + pG T Y y 



QH 

X rQH - kY 



(8) 

It is worth noting that typically Ty/G <C 1, which means 
that the viscous force is usually much weaker than the 
gyrotropic one H, Therefore, a good starting point 
would be the frictionless limit Fy = 0. In that case 
the vortex moves along the lines of constant potential 
U(X, Y) = const. From that one can deduce a crossing 
time T = irj (■jfioH) that is remarkably ^insensitive to the 
detailed structure of the domain wall [13], as indeed ob- 
served experimentally [l8[. However, the viscous loss of 
energy is a crucial factor determining the average veloc- 
ity of a domain wall: any drift reflects the dissipation of 
the Zeeman energy —QHX; in the frictionless limit the 
wall exhibits no drift at all. Thus one must include the 
effects of viscous friction to evaluate the drift velocity. 
A general solution of the equations of motion {5} reads 



X 
Y 



Y{pG-T X Y)/Txx = Vt 



const, (9) 
(10) 



where n = (G 2 + detT)/{kT X x) 
-(p-X9)GQH/(kT X x), and g 
Two distinct regimes are found. 



w G 2 /(kT XX ), = 
- (rV X x - x^xy)/G. 
At low applied field, 
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FIG. 2: The drift velocity Vd of the domain wall as a function 
of the applied field H for a permalloy strip of width w = 200 
nm and thickness t — 20 nm. Dashed vertical lines mark the 
critical fields H c - and H c + . Symbols are results of numerical 
simulations with in-plane mesh sizes as shown. 



the equilibrium position of the vortex is inside the strip. 
After a relaxation period of duration t\ ~ G 2 /(kTxx) 
the wall reaches a state of steady drift with X = V = 
HlfH (/j-lf = Q/^xx is the mobility in low fields), and 
Y = Yoo ~ —pGV/k. Note that in the absence of the 
gyrotropic force, the relaxation time would have been 
much shorter, Tyy jk. The gyrotropic effect is apparently 
one of the reasons why the mode £i = Y is particularly 
soft. 

Above a critical field the restoring potential fails to 
prevent the vortex from reaching the edge, where it 
merges with the antihalfvortex. Our numerical experi- 
ments (see below) indicate that the vortex is immediately 
re-emitted with the same chirality x an d opposite polar- 
ization p and starts to move towards the opposite edge 
(Fig. 1, bottom). The critical fields are slightly different 
for p — +x and p = H c ± = H c0 /(l =p g), where 

H c o = ^hFkw/2G and j « 1. In the narrow interval 
i/ c _ < H < H c+ the vortex reaches a steady state for 
p = +x but not for p = — x- As one might expect, the 
breakdown of steady motion coincides with the softening 
of the first mode: at H = H c $ the crossing time T = 2t\ . 

Above H c+ the vortex crosses the strip regardless of 
its polarization, and an oscillatory regime sets in. For 
the drift velocity 14 we find 



V d = AiLF# 



214(1 + detr/G 2 



(11) 



atanh(# c+ /ff) + a,tanh(H c - / H) 

At first, the drift velocity drops precipitously (Fig. [5]), 
changing its order of magnitude from 0(a^ 1 ) to 0(a). 
In higher fields the velocity once again becomes propor- 
tional to H, albeit with a smaller mobility /iHp: 



A'HF 



(r 2 T xx 



2rxT X Y + r y y)r xx 



G 2 



< I. 



(12) 



For a quantitative analysis [3| we turn to the model of 
a vortex domain wall of Youk et al. ■ The composite 



wall consists of three 90° Neel walls comprising the anti- 
halfvortices and a vortex that can slide along the central 
Neel wall (Fig. [T]). We used saturation magnetization 
M = 8.6 x 10 5 A m" 1 , Gilbert damping a = 1CT 2 , and 
exchange constant A — 1.3 x 10~ n J m _1 , yielding the 
exchange length A = 'A/ 'noM 2 =3.8 nm. 

The damping coefficients ([4]) are determined mostly 
by areas with a large magnetization gradient Vm, i.e. 
from the three Neel walls whose width is of order the ex- 
change length A, which gives Tij ~ aJtw/X. The values 
of damping coefficients are as follows (l9j : 

T X x = 0.044G, T XY = 0.031xG, T YY = 0.049G. (13) 

The stiffness constant k of the restoring potential could 
not be calculated accurately because two of its main con- 
tributions, a positive magnetostatic term and a negative 
term due to Neel-wall tension, nearly cancel out. This 
is not surprising given the proximity to a region where 
the vortex wall is unstable [5[ . Instead, we extracted the 
relaxation time t\ directly from the numerics (see be- 
low) by fitting Y(t) to Eq. (10|). We obtained n in the 
range from 8.5 to 9 ns for fields from 4 to 60 Oe with 14o 
scaling linearly with H. In calculating the critical veloc- 
ity 14 = kw/(2G), we replaced w with an effective strip 
width Wes = w — 2R, where R is a short-range cutoff 
due to the finite size of a vortex core [12]. From vortex 
trajectories observed numerically (top panel of Fig. [3]) we 
estimate R ks 10 nm. 

To compare our theory with experimental results, we 
have computed the low and high-field mobilities using 
standard material parameters for permalloy (see meth- 
ods) for a strip of w = 600 nm and t = 20 nm employed 
in the experiment of Beach et al. [jo| . While the cal- 
culated low-field mobility /z^p = 29 m s _1 Oe~ 1 agrees 
reasonably well with the experimental result = 
25 m s _1 Oe _1 , our estimate of the high-field mobility 
^hf = 0-61 m s _1 Oe _1 is markedly lower than the ob- 
served value fi^p — 2.5 m s _1 Oe _1 . 

To understand the discrepancy between theory and 
experiment at high fields, we compared the theoretical 
curve Vd(H) against numerically simulated motion of 
a vortex domain wall in a permalloy strip with width 
w = 200 nm and thickness t — 20 nm. Numerical simu- 
lations were performed using the package oommf [2l|. We 
used the same material parameters as mentioned above. 
Cell sizes were 2 nm x 2 nm x 20 nm for most runs and 
5 nm x 5 nm x 20 nm in a few others. The strip length 
was L = 4 /zm or more. Care was taken to minimize the 
influence of a stray magnetic field created by magnetic 
charges at the ends of the strip. 

The drift velocity 14 computed within the two-mode 
approximation agrees reasonably well with simulation re- 
sults both below and above the breakdown field H c+ = 
9.5 Oe up to a field of H 2 « 35 Oe (Fig. EJ. However, 
above Hi the numerically observed drift velocity begins 
to increase in disagreement with the theory. The failure 
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FIG. 3: Top: The transverse vortex coordinate Y(t) for sev- 
eral values of the applied field H. Deviations from the ex- 
pected behavior (|10|l in weak fields are due to stray field from 
the strip ends. Bottom: The width of the wall A(t). Curves 
for different fields are shifted vertically by 150 nm for clarity. 
The initial width in all cases was A(0) = 190 nm. 



of the two-mode approximation around Hi was traced to 
the softening of another mode seen as fast oscillations of 
the width of the domain wall (the width was measured 
as the difference in x-coordinates of the antihalfvortices, 
top panel in Fig. [3]). The new mode is excited at the 
beginning of each cycle and relaxes to an equilibrium on 
the time scale T2 ~ 2.5 ns. In a field of H = 24 Oe this 
mode decays well before the end of the cycle (T = 7.4 
ns, see the bottom panel of Fig. [2]). It is responsible 
for a small fraction, Ofo/T), of the net energy loss and 
thus can be neglected. At H = 48 Oe (T = 3.7 ns) the 
new mode stays active all the time and therefore cannot 
be ignored. In accordance with this, the numerical data 
begin to deviate from our two-mode model (|11| around 
H-2 = 35 Oe. The new mode is related to the incipient 
emission of an antivortex by one of the edge defects. A 
similar mechanism may be at work in wider strips used 
by Beach et al. 

The framework presented here is sufficiently simple and 
flexible to include additional modes and the effects of spin 
torque. It can also handle other scenarios observed in 



numerical simulations wherein the absorbed vortex is re- 
emitted with the opposite chirality [ItJ or not re-emitted 
at all Q or the vortex core flips while the vortex is still 
in the bulk 0, HI). Antivortex walls @, E3, H3 can be 
handled in a similar way, provided one develops a simi- 
larly detailed model to compute the energy and damping 
coefficients. 
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